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ABSTRACT 



Imaging spectroscopy of the solar atmosphere using the Na i D[ line yields marked asymmetry between the blue and red line wings: 
sampling a quiet-Sun area in the blue wing displays reversed granulation, whereas sampling in the red wing displays normal granula- 
tion. The Mg i b2 line of comparable strength does not show this asymmetry, nor does the stronger Ca n 8542 A line. We demonstrate 
the phenomenon with near-simultaneous spectral images in Nai Di, Mgibi, and Can 8542 A from the Swedish 1-m Solar Telescope. 
We then explain it with line-formation insights from classical ID modeling and with a 3D magnetohydrodynamical simulation com- 
bined with NLTE spectral line synthesis that permits detailed comparison with the observations in a common format. The cause of 
the imaging asymmetry is the combination of correlations between intensity and Dopplershift modulation in granular overshoot and 
the sensitivity to these of the steep profile flanks of the Na i Dj line. The Mg i b^ line has similar core formation but much wider wings 
due to larger opacity buildup and damping in the photosphere. Both lines obtain marked core asymmetry from photospheric shocks 
in or near strong magnetic concentrations, less from higher-up internetwork shocks that produce similar asymmetry in the spatially 
averaged Can 8542 A profile. 
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1. Introduction 

Three technological advances currently improve optical stud- 
ies of the solar atmosphere. First, adaptive optics and numer- 
ical postprocessing deliver angular resolution in groundbased 
imaging well beyond the atmospheric Fried limit (rarely larger 
than 10 cm equivalent aperture even at the best sites). Second, 
multi-line imaging spectroscopy leaped forward with fast-tuning 
narrow-band Fabry-Perot interferometers at the major solar tele- 
scopes. Third, numerical MHD simulations have evolved into 
productive analysis tools of small-scale near-surface magneto- 
convection and are becoming reliable emulators of fine structure 
higher up in the solar atmosphere. 

In this paper we combine these advances in a comparative 
study of the formation of three of the strongest lines accessible 
to the current generation of Fabry-Perot instruments, which all 
observe in the red part of the spectrum. The three lines are the 
b 2 line of Mgi at A = 5172.698 A (henceforth "Mg line"), the 
Di line of Nai at A = 5895.940A (henceforth "Na line"), and 
the infrared-triplet line of Can at A = 8542. 144 A (henceforth 
"Ca line"). The wa velengths are for standard air and taken from 
iMooreet al.l ( 19661) 

These three lines are often taken to be diagnostics of the solar 
chromosphere, together with Ha which gave the chromosphere 
its name through its off-limb spectral dominance (Lockver 1868 
see iRuttenl 2010). On the solar disk Ha displays dense canopies 
of fibrils that cover active regions, network, and even quiet-Sun 
internetwork cell interiors except in the very quietest regions. 
The detailed comparison between Ha- and the Ca line by Cauzzi 
et al. (2009) shows that the latter line displays these chromo- 



spheric fibrils similarly near quiet-Sun network but with smaller 
extent, whereas in quiet internetwork cell centers it samples 
the domain of generally cool but shockwave-ridden gas below 
the fib rilar canopies that was termed "clapotisphere" by IRuttenl 
(1995) We use this name here also, reserving "chromosphere" 
for the fibril canopies seen in Ha. The Na core is mostly formed 
within the wave-domin ated clapotispheric domain, as shown by 
iLeenaarts et al.l (20101) We show below that the same holds for 
the Mg core. 

The present study i s a companion to the observational Ha- 
Ca line comparison of lCauzzi et ail (20091) by now comparing 
the Na and Mg lines to the Ca line. It is a direct extension 
of the observation-simulation comparisons by Leenaarts et al.l 
(20091) for the Ca line and by ILeenaarts et al.l (20101) for the 
Na line, by now comparing these lines with another and with 
the Mg line. However, in this study the emphasis does not lie 
on the chromospheric or clapotispheric response of the cores 
of these lines but rather on their inner-wing formation in lay- 
ers only 100-200 km kilometer above the surface, where the 
phenomen on of "reversed granulation" dominates the observed 
scene (e.g., Leighton et al . 1962, Evans 1964, E vans & Cat alano 
19721; ICan field & Mehltrettej 19731 lAltrock & Musmanl 197 1 
Nordlund 1984b 



Suemoto et al. 1987 1990, Rutten et al. 200 



Leenaarts & Wedemever-Bohm 2005 , Janssen & Cauzzi 2006; 
ICheung et al. 2007 ). The topic of this paper is not its physical 
nature, but how how this domain appears in the three lines. 

This study started with the data shown in Figure Q] from 
the CRISP imaging spectrometer at the Swedish 1-m Solar 
Telescope (SST, La Palma). They sample the Mg and Na lines 
from a small active region. We first inspected the line-center im- 
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Fig. 1. Profile-sampling image comparisons between the Mg line (upper row) and the Na line (lower row). These images are 
subfield cutouts from a dual-line SST/CRISP profile scan taken on August 23, 2010. They sample a small active region with pores 
in the upper part and a quiet area in the lower part of the subfield. Each image is byte-scaled independently and is labeled with its 
wavelength separation AA from line center in mA. Both sequences sample the line profile from the outer blue wing through the core 
to the outer red wing. The AA values in the first two columns were selected for scene similarity between the two lines. The next two 
columns were also selected for scene similarity but have equal AA values for the two lines. The final three columns mirror the AA 
values of the first three. All images show rich detail; we invite the reader to magnify them with a pdf viewer. This paper explains 
the imaging asymmetry between the two lines in the red-wing panels where the Mg line shows reversed granulation as in its blue 
wing, whereas the Na line instead shows normal granulation as in the last panel. 



ages in the fourth column. They display bright fine-scale fibrilar 
structure near some pores, next to dark internetwork elsewhere. 
Thus, the Mg and Na cores do display some chromospheric re- 
sponse, i.e., can be opaque in fibrils. These line-center images 
are remarkably similar for the two lines. 

We then compared the images taken in the blue wings of the 
two lines and quickly found that for every scene observed in the 
one line, a nearly identical scene can be found in the other. One 
only has to search further out into the much stronger Mg wing. 
We concluded that these wings sample the same solar layers sim- 
ilarly, and that there is no point in observing them both. 

However, when we then turned to the red wings we were sur- 
prised by the striking difference illustrated in Figure Q] in the red 
Na wing the reversed granulation drops out. The Mg line shows 
reversed granulation in its red wing that is similar to that in its 
blue wing, but in the Na red wing the line-center scene of chro- 
mospheric activity is quickly replaced by granulation, without 
an intermediate stage exhibiting reversed granulation. Along the 
red Na wing one seems to suddenly drop from chromospheric 
and clapotispheric sampling at line center to photospheric gran- 
ulation sampling further out. 

These inspections were made at the SST while T.D. Tarbell 
(LMSAL, Palo Alto) and P. Siitterlin (ISP, Stockholm) 
were present. The first commented that he had noted this 
reversed/non-reversed granulation asymmetry between the Na 
line wings in Hinode/SOT images and had wondered about its 
cause. The second suggest ed granular brightnes s-Dopplershift 
correlation as in Figure 6 of lSiitterlin et al.l (2001l). 

Their suggestions were followed up in an observing cam- 
paign at the SST during October 2010, targeting truly quiet Sun 
to better view reversed granulation and adding the Ca line be- 
cause it must sample the same scenes somewhere along its ex- 
tended wings. We analyze the best resulting data here, com- 
paring them to modeling predictions from a three-dimen sional 
MHD simulation snapshot following Leenaarts et al. (2009, 
201011. 

Below, we describe the observations in Section [2] the sim- 
ulation and the spectral synthesis in Section [3] and the results 



in Section [4] Since the latter consist of profile-sampling image 
cubes for the three lines for both observations and simulations, 
we use a common format permitting direct comparison to inspect 
the observed and simulated cubes in a parallel descriptive char- 
acterization of what they show. We then start the interpretation 
in Section |5] by reviewing what one expects from mean-profile 
inspection and from NLTE line synthesis with a standard one- 
dimensional hydrostatic solar-atmosphere model. We use these 
insights to interpret the results and explain most observed fea- 
tures. We add Hinode images as indicator of center-limb behav- 
ior, and end the paper with discussion in Section Q and conclu- 
sions in Section[8] 

2. Observations 

2. 1 . SST/CRISP observations 

The C Risp Imaging SpectroPolarimeter (CRISP. [Scharmer et al. 
2008 ) at the Swedish 1-m Solar Telescope (SST. IScharmer et al. 
2003a) contains a dual Fabry-Perot tunable filter system. It is 
designed for diffraction-limited narrow-band imaging with fast 
passband tuning sampling multiple spectral lines at rapid ca- 
dence in the wavelength range 5100-8600 A. 

The observing campaign, which was educationally oriented 
involving Oslo students on-site, suffered from bad observing 
conditions. However, on the last day, October 30, 2010, there 
was a 40-min period of moderate-quality seeing permitting ob- 
servation of a region close to disk center (Stonyhurst helio- 
graphic coordinates (N4, W0) degrees) that was selected for its 
absence of any rosettes or other chromospheric activity in Ha. 

The observing sequence sampled four spectral lines: the Mg, 
Na and Ca lines discussed here and also Fei 6302 A at AA = 
-0.048 A from line center to obtain magnetograms. The Mg line 
was sampled at 41 wavelengths between AA = +0.892 A from 
line center with a passband with FWHM = 32 mA, the Na line at 
41 wavelengths between AA = ±1.715 A with passband 57 mA, 
the Ca line at 47 wavelengths spanning AA = +2.717 A with 
passband 111 mA. The sampling had regular spacing across each 
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core but was sparser in the outer wings. The scan durations were 
about 10 s per line. At every wavelength step 8 exposures were 
acquired. Since the cameras operated at 35 frames/s speed and 
the prefilter change time ranged from 0.27 s to 0.52 s, the total 
duration for one complete scan through all 4 lines was 29.7 s. 
We selected the best scan for this analysis. It was taken from 
09:09:03 to 09:09:33 UT. 

The spatial resolution of these observations was improved 
beyond the seeing Fried li mit by using the S ST ada ptive op- 
tics system described by I S charmer et al . (2003b) and by 
subsequent postprocessing with the Multi-Obj ect Multi-Frame 
Blind Deconvolution (MOMFBD) technique of lvan Noort et alJ 
(2005) M ore details on such image re d uction are given in, e.g., 
Ivan Noort & Rouppe van der Voortl (2008) and Rouppe van 
der Voort et al. (2009) The field of view after this process- 
ing and co-alignment is 58 x 56 arcsec 2 , with 0.0592 arcsec/px 
or 42.6 km/px scale. However, the figures below show only a 
smaller area which equals the horizontal simulation extent. 

2.2. Conversion to brightness temperature 



The brightness temperature 7b is the formal temperature that 
yields fi,i(7b) = 7° bs , i.e., that reproduces the observed intensity 
when entered into the Planck function. It equals the gas tem- 
perature Tgas at the height where — 1 along a radial line of 
sight when LTE and the Eddington-Barbier approximation are 
valid. LTE does not hold for the cores of our lines, but we never- 
theless convert intensities into brightness temperatures to enable 
comparisons that are not affected by the differences in Planck 
function sensitivity to temperature at different wavelengths. 

The inverse Planck operation requires rescaling of the ob- 
served intensities into absolute units. We did this by compari- 
son to the solar disk-center spectrum atlas obtained by Brault 
and Testerman with the F ourier Transform Spect rometer (FTS) 
at Kitt Peak, c alibrated bv lNeckel & Labs! (1984) and spread by 
INeckell dl999l) . We smeared the atlas profiles of the three lines 
with the corresponding CRISP spectral transmission functions, 
i.e., the product of the measured prefilter transmission and the 
theoretical transmission of each of the two etalons at the line 
wavelengths. The resulting functions have a prominent symmet- 
rical sidelobe pattern, with the first pair reaching 0.001 relative 
amplitude at AA = +2.2 A for the Na line. Convolving the atlas 
profiles with these makes the line cores rise substantially. For 
each line, we rescaled the data units of the spatially-averaged 
observed profile to the absolute intensities of the smeared at- 
las profile by determining their ratio for the wavelength with 
the highest outer-wing intensity in the observed profile. The ob- 
served and the smeared atlas cores show satisfactory agreement 
(last column of Figure|9j, confirming that the transmission func- 
tions give an adequate description of the spectral smearing. 
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Fig. 2. Intensity contrast comparisons for the three lines 
between observations (solid), simulation (dashed), and the 
smeared simulation (dotted). Each curve shows the normalized 
occurrence distribution of pixels with the given contrast (inten- 
sity divided by the spatial mean) for the sum of the blue and red 
outermost-wing image pairs. 




10 15 

x [arcsecl 



20 



5 H 



10 



15 



x [Mm] 



Fig. 3. Unsigned magnetic field distributions. Top: observed 
Stokes- V amplitude at AA = -0.048 A in the Fei6302A line 
for the subfield. Only a few small magnetic concentrations have 
signal above the noise. The two major ones near x — 10 have 
opposite polarity. Bottom: field amplitude in the simulation at 
h = km, corresponding to spatially-averaged T5 = 1 in the con- 
tin uum at A — 5000 A. The field is bipolar, as shown in Figure 5 
of iLeenaarts et al.l (20101) Throughout this paper image scales 
are in arcsec, simulation scales in Mm. The two areas have equal 
solar-surface extent. 



2.3. Hi node observations 

The Hinode satellite (iKosud etalj 2007) was programmed to ob- 
serve at several center-to-limb positions along the central merid- 
ian on December 4, 20 1 between 1 4: 1 4: 3 3 and 1 5 : 3 3 :4 1 UT us- 
ing the Narrowband Filte r Imager (NFI) fed by the 0.5 m Solar 
Optical Telescope (SOT, iTsuneta et al.l 20081; ISuematsu et alJ 
20081). Successive scans sampling the spectrum in 80 mA steps 
from -480 to +480 mA from line center were made through the 
Mg line and the Na line, with bandwidths 74 and 84 mA, re- 
spectively. One line scan took approximately 22 s. The prefilter 
changes and context images taken with the Broadband Filter 



Imager added 10 s delay between successive line scans. The data 
were processed with the standard SolarSoft fg_prep pipeline. We 
selected small subfields of the size of the simulation extent from 
the pointings at p — 1 .0, 0.5, and 0.35, where p = cos 6 with 6 the 
viewing angle, for display in Section|6] 

3. Simulations 

3.1. Simulation snapshot 

We employ the same snapshot from a 3D radia tion-MHD sim- 
ulation performed with the Oslo Stagger Code (lHansteen et al.l 
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Fig. 4. Observed brightness-temperature im- 
ages for selected wavelengths sampling the 
profiles of the Mg line (first column), Na line 
(second column), and Ca line (third column). 
The area is a small cutout of the center of the 
full field of view that corresponds in size to the 
horizontal extent of the numerical simulation. 
The identifier in each panel specifies the line, 
the CRISP passband offset from line center 
in km/s with redshift positive, and the mean 
brightness temperature over the full field of 
view. The grayscale range of each panel is 
clipped at AT b = ±400 K relative to this mean, 
so that each panel covers the same brightness 
temperature range in magnitude. The passband 
offsets select equal FALC formation height 
across the first and last two rows (optical depth 
unity near h = 65 and 150 km, respectively) 
and near-equal wavelength offset AA from the 
nominal line center wavelength in Dopplershift 
units across the middle three rows (about 
AA = -3, 0, and +3 km/s, respectively). The 
selections are symmetrical with respect to line 
center (central row). 
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2007) that w as used inlLeenaarts et alj (20091) for Ca line syn- 
thesis and in lLeenaarts et al.l (20101) for Na line synthesis. This 
simulation snapshot contains significantly more magnetic field 
than the observations (Figure[3]), but we prefer to re-use the snap- 
shot and spectral line synthesis nevertheless because otherwise 
we would have to re-do and re-dis play much a nalysis and diag- 
nosis presented in Leenaarts et al. (2009 2010 1 that we can now 
simply refer too. 

The (x,y, z) snapshot grid measures 256 xl 28 x 160 sampling 
points, corresponding to a physical size of 16.6x8.3 x 15.5 Mm 3 
with the height ranging from the upper convection zone to the 
corona. The grid spacing is 65 km/px in the horizontal x and 
y directions; in the vertical z direction it varies from 32 km at 
the bottom to 440 km at the top. The snapshot contains bipolar 
magnetism with a mean field strength of 150 gauss in the photo- 
sphere. The vertical field magnitude at t he bottom of the photo- 
sphere is shown in Figure [3] Figure 5 of lLeenaarts et al.l (20101) 
displays various physical quantities, including field strength, at 
different heights in the simulation. 



pier profile synthesis of the Mg line was added by using MULTI 
and assuming each vertical column of the simulation snapshot to 
describe the stratifications of a plane-parallel ID atmosphere. 

Because 3D radiative transfer computation is very com- 
putationally demandin g, the model atom for Na described in 
iLeenaarts et al.l (2010) contains only the levels governing the 
line, the continuum, and an intermediate level with large colli- 
sion bound-free coupling to ensure the "photon sucti on" popula- 
tion flow that tends to populate the Na i ground state dBruls et al.l 
1992). The model atom for Ca was a simple 5 -level plus con- 
tinuum one which is adequate to describe the formation of Ca n 
H&K and the C an infrared triplet including our Ca line (e.g., 
lUitenbroe k 1989). The model atom for Mg used he re was a 13- 
level plus continuum model that was also used by iFleck et al.l 
20 101 and is based on th e 65 -level plus continuum model atom 
of ICarlsson et al. ( 19921) The level reduction used the method 



of lBard & Carlssonl (2008) Photoionization was evaluated with 
line-blocking accounted for by opacity samplin g in 9000 wave- 
length points as described hy Collet et al. 20051 



3.2. Radiative transfer 

For the Ca line and the Na li ne we also r e-use the line synthesis 
described in Leenaarts et al. (2009 2010 ), respectively, to which 
we refer for further detail. These radiative transfer computa- 
tions were perfo rm ed in 3D geometry using the c ode Multi3D of 
iBotne n ( 1997) and Leenaarts & Ca rlsson (2009) which is based 
on the ID code MULTI of lCarlssonl ( 19861) For this paper, sim- 



3.3. Spectral and spatial smearing 

In order to mimic the observations we first smeared the simula- 
tion intensities in the spectral domain with the CRISP transmis- 
sion functions discussed in Section 12.21 Their effect is large in 
the line cores and small in the outer wings (Figure |9}. 

Subsequently, we smeared the simulation spatially with an 
estimated SST point spread function (PSF). Its determination is 
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Fig. 5. Simulated brightness-temperature im- 
ages for selected monochromatic wavelengths 
sampling the profiles of the Mg line (first 
column), Na line (second column), and Ca line 
(third column). The area equals the observed 
subfield in Figure [4] in size. The format, 
labels and grayscales correspond to those of 
Figure [4] i.e., each panel is again clipped at 
AT/j = ±400 K relative to the spatially aver- 
aged mean. The wavelength selection, again 
specified in Dopplershift units, corresponds 
to the selection in Figure [4] i.e., near-equal 
FALC formation height of about 65 km and 
150 km along the outer rows and near-equal 
wavelength offset AA from the nominal line 
center in Dopplershift units across the middle 
three rows, of which the center one samples 
nominal line center. 
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less straig htforward. Even the shape of the typical PSF remains 
in debate. Levvl (1197 11) 1 adop ted t he sum of two Ga ussians, as 
did iDeubner & Mattid ( 19751) but iNordlundl ( 1984a|) advocated 
the sum of two Lorentzians instead. Collad os & Vazquezh 19871) 
tested both formalisms and concluded that the choice strongly 
affects granulatio n contrast values that r esult from restoration 
by deconvolution. ISchiissler et al.l (20031) used the sum of the 
SST's theoretical Airy function and a Lorentzian to describe 
the atmosphe ric degradation, while the comparable paper by 
IShelvag et al.l (2004) smeared the same simulation with the sum 
of the DOT's theoretical Airy func tion and a two-parameter win g 
function following Equation 6 of ICollados & Vazquez! (19871) 
which is not a Lorentzian but describes the straylight across a 
knife-edge such as the lu nar limb when the PSF is Lorentzian. 
IWedemev er-Bohm ( 2008) and Wedemeyer-Bohm & Rouppe van 
der Voort (2009) defined PSFs for seeing-free imaging with 
Hinode/SOT as the convolution of the theoretical Airy func- 
tion and a Voigt function, best-fitting the Voigt parameters with 
Mercury transit and eclipse data. 

In our case both seeing and scattering cause angular smear- 
ing, we have no eclipse or transit data to measure their ef- 
fects, and we must cope with the wavelength variation be- 
tween our lines. We therefore do not aim to obtain a reliable 
PSF for data deconvolution, but simply apply somewhat real- 
istic smearing to the simulation as illustration of its likely ef- 
fect. We chose to use a single Voigt profile, representing the 
convolution of the Gaussian approximation to the Airy function 
A(r) w l/( V2~7r w) exp(-r 2 /2w 2 ) with extended Lorentzian tails 



L(r) = y/(4 7r 2 )/[(y/47r) 2 + r 2 ] describing scattering, and assume 
circular symmetry with radius r as variable. In principle, the Airy 
width is w — 0.42 A/D, but since the seeing during this observa- 
tion was only mediocre (in terms of the superb seeing that the La 
Palma site offers at its best), the MOMFBD restoration probably 
did not reach the full-aperture resolution. We therefore reduce 
the aperture D entering the Airy component to an "effective" 
one with a fractional factor / and scale the res ulting w oc A~ 1 ^ 5 
as given by atmospheric turbulence theory (see iRoddierl 1981 ). 
The Voigt parameter a = y/(4n V2 w) is scaled with A s . We thus 
have three fit parameters: /, a, and s. 

As fit constraints we used the resulting contrast distribution 
and the visual appearance of the granulation at the extremes of 
our common wavelength ranges, i.e., the deepest-formed wave- 
lengths that permit observation-simulation comparison. These 
differ per line. We use these extremes because the simulation 
is probably most realistic at the granulation level (indeed, most 
recent solar PSF estimates aim to show that granulation sim- 
ulations are more reliable than PSF-challenged observations). 
We sum these outer wing images as a first-order compensation 
for Dopplershifts. (It would be better to use parallel nearby- 
continuum observations but the CRISP prefilters do not permit 
taking such.) We found that the combination of / = 0.5 and 
a = 0.2 at the wavelength of the Mg line with s = —2 pro- 
duces both the acceptable contrast fits shown in Figure [2] and 
similar appearance of the observed and smeared granulation (not 
shown, but compare the outer rows of Figures|4]and|6}. The value 
/ = 0.5 seems realistic for non-superb seeing; a = 0.2 adds ex- 
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Fig. 6. The same simulated brightness- 
temperature images as in Figure [5] but after 
spectral smearing with the CRISP transmission 
profiles described in Section 12.21 and spatial 
smearing with the point spread function esti- 
mates obtained in Section[3~3l 
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tended wings below the 1% level, and s = -2 describes scattering 
between the Rayleigh and Mie limits. The differences in distri- 
bution width between panels in Figure [2] primarily reflect dif- 
ferent Planck function sensitivities to temperature. The dip and 
enhanced righthand tail, most clearly present for the Mg line, 
are due to the appreciably larger amount of magnetism in the 
simulation. 



4. Results 

In this section we present displays from the profile-sampling im- 
age cubes for the three lines from the observations and from 
the simulation, likewise and in a common format. These fig- 
ures serve to supply a comparative empirical inventory of what 
quiet-Sun fine structures the three lines display, in parallel for 
observation and simulation, while postponing all interpretation 
to Section [5] We again invite the reader to magnify the figures 
per pdf viewer for better appreciation of their detail. 



4.1. Magnetograms 

Figure[3]compares the magnetic field in the subfield of the obser- 
vations and in the simulation. The observed magnetogram indi- 
cates apparent longitudinal flux density, the simulation panel dis- 
plays intrinsic vertical field strength. The observed subfield has 
only a few very small magnetic concentrations; there are some 
larger ones elsewhere in the full field of view. The field in the 
simulation is shown quantitatively in Figure 5 of iLeenaarts et al.l 



(2010). displaying the polarity distribution in the first panel and 
the lateral expansion of the magnetic concentrations with height 
in subsequent panels. The larger white patches in the simula- 
tion panel represent kilogauss concentrations. The few observed 
white patches are likely to also represent kilogauss fields. There 
is much more strong field in the simulation than in the observa- 
tion, but it is nevertheless a relatively weak-field simulation; the 
field does not upset granular convection as is the case in strong 
network or plage. 

4.2. Profile-sampling images 

Figure|4]shows a selection of the observed line-sampling images, 
cut down to a small subfield at the center of the SST field of 
view of which the 23.1 x 1 1.5 arcsec 2 area equals the horizontal 
simulation extent. All panels have dynamic range AT/, = +400 K 
around the mean value over the full field of view. The conversion 
into brightness temperature and the common greyscale range 
make all panels intercomparable. Similarly, the wavelength sep- 
arations from line center AA are specified in Dopplershift units 
(km/s) to supply wavelength-independent comparison. 

The AA selections are specified in each panel and are sym- 
metric along the columns around the line-center panels in the 
central row. They have been carefully selected to obtain similar 
atmospheric sampling in the three lines. For the first and second 
rows (and reversely the last and next-to-last rows) the AA val- 
ues are those that produce equal r,i = 1 formation height in the 
FALC modeling described in Section IBT21 below, at h = 65 km 
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Fig. 7. Upper assembly: observed brightness- 
temperature Dopplergrams for the three lines 
in the format of Figure [4] Each panel displays 
the normalized Dopplergram subtraction 
D = (T™' - T* lw )/(T^ + r b bk,c ), with bright 
implying blueshift (updraft) for spectral ab- 
sorption features and redshift (downdraft) for 
spectral emission features. The grayscale is 
defined so that D = is represented by the 
same gray in all panels while the ranges are 
adapted per panel for largest contrast. The 
separations from line center AA, specified in 
km/s in each panel, are the same as in Figure|4] 
Lower assembly: corresponding results from 
the smeared simulation. The AA values are the 
same as in Figures[5]and[6] 
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and h — 150 km, respectively. The first height is the lowest at 
which we can do such three-line comparison, being limited by 
the outermost CRISP tunings in the extended wings of the Ca 
line. This height samples the top of the normal granulation. The 
second height was selected to sample reversed granulation. At 
both heights LTE is a good assumption (Section l5T2l i; indeed, the 
mean 7j values specified in each panel are similar along these 
rows and between opposite-wing pairs. Closer to line center we 
instead selected similar AA values along rows, since these im- 
ages may primarily sample clapotispheric structures in this very 
quiet area. Cloud-like line formation with specific cloud Doppler 
velocities should then produce similar scenes for lines with sim- 
ilar cloud opacities and source function variations. 

The granulation appears closely similar along the first row, 
but not so along the bottom row where the Na line displays the 
largest granulation contrast, the Ca line the smallest. The re- 
versed granulation also appears closely similar in all three lines 
along the second row, but not so in the next-to-last row where 
the Na line roughly shows the same granulation as in the bot- 
tom panel, whereas the Ca line displays the closest similarity to 
its blue-wing counterpart in the second row. These differences 
confirm the imaging asymmetries in Figure Q] 

The line core scenes in the central three rows are similar for 
the Mg and Na lines, but differ for the Ca line. For all three lines 
they differ much on the opposite sides of line center, with oppo- 
site grayscale patterns especially for Mg and Na. These patterns 
are patchy in morphology and appear fuzziest in the Na line. In 
the Ca line small roundish bright features appear, clearest on the 
blue side but not on the red side. 



The few magnetic concentrations in this very quiet area 
(Figure [3} appear as bright points, clearest in the second row 
(all three lines) and next-to-last row (Mg and Ca). 

Figure [5] presents the companion display from the simu- 
lation. The AA selections are the closest available to those in 
Figure |4] In all panels the fine structure has appreciably larger 
brightness-temperature contrast than in Figure particularly in 
the central three line-core rows. The granulation has largest con- 
trast along the top and bottom rows in the Mg line, smallest in 
the Ca line. The reversed granulation is most clear for Na in 
the second row, least clear for Na in the next-to-last row. In the 
latter it appears present as well, but all the bright patches and 
lanes in this panel represent strong magnetic field (compare with 
Figure|3]and the bottom Na panel). 

In the three line-core rows of Figure |5] the Mg line dis- 
plays dark intergranular lanes in all three rows, most strikingly 
at AA = +3.6 km/s. They are covered by higher-lying structures 
in the other two lines. These are grayish with much fine struc- 
ture at Na line center but become black-and-white clipped at 
the 400 K display limits at Ca line center. They break up into 
smaller patches in the near-center panels, many opposite in black 
or white at the two sides. 

Figure[6]repeats Figure |5]after smearing. The simulation pix- 
els (65 km) are actually coarser than the SST pixels (42.6 km 
on the observing date), but comparison of Figures |4] and [5] sug- 
gests that the actual resolution of the observations is consider- 
ably worse. We therefore apply spectral smearing with CRISP's 
transmission profiles and spatial smearing with the estimated 
PSFs from Section [3~3l to illustrate their effect. The central three 
rows of Figure [6] are most affected, especially with respect to 
the extended patches in the Na line. At line center, these loose 
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Fig. 8. Scatter diagrams of the core asymmetry 
ratio V/R = T™ w /T£* against field strength. 
Upper two rows: observations, using the full 
field of view. The abscissa is the magnitude of 
the apparent longitudinal flux density measured 
from the Fe 1 6302 A magnetogram, in arbitrary 
units and after 9-pixel spatial smoothing to 
mimic field spreading with height. Bottom row: 
smeared simulation. The abscissa is the mag- 
nitude of the vector field at height z = 500 km. 
Left: Mg line. Center: Na line. Right: Ca line. 
The sampling wavelength in km/s is specified 
in each panel. The V/R values are multiplied 
by 1000 for plotting reasons. Density contours 
replace individual pixels to avoid plot satu- 
ration. The solid curves are the distributions 
over V/R and field strength, respectively. 
The dashed curves show the first moments 
per bin per axis. They line up for high cor- 
relation and are perpendicular at no correlation. 



the fine structure seen in Figure |5] and become smoother grey- 
ish patches, appreciably darker than the bright magnetic concen- 
trations and with similar appearance as the extended patches in 
the corresponding panel of Figure |4] Dark intergranular lanes 
remain visible in the Mg core. The Na-line imaging asymme- 
try of reversed granulation in the second and next-to-last rows 
remains, but without showing dark lanes outside the field-filled 
ones in the red wing as observed in Figure |4] 

4.3. Dopplergrams 

The various asymmetries above between the blue and red line 
wings obviously suggest a role of Dopplershifts. Therefore, 
Figure|7]shows Dopplergrams, in the upper half for the observa- 
tions and in the lower half for the smeared simulation. The sign 
of the Dopplergram definition D = (7^ ed - 7* lue )/(T£ :d + r£ lue ) is 
chosen to make the granular Dopplershift pattern appear as gran- 
ulation, i.e., with bright granular updrafts and dark intergranular 
downdrafts. The D ranges are optimized per panel, but D — is 
the same dark gray in each panel. The AA selections correspond 
to those in Figures |4]-|6] The first rows again sample the top of 
the normal granulation, the second rows the reversed granula- 
tion, the third rows the inner line cores. 

In the observed Dopplergrams the pattern quality along the 
first row shows that the Na line has the best granular Doppler 
sensitivity, Ca the worst. Most structure in the Ca Dopplergram 
is noise or MOMFBD artifact that is much exaggerated by the 
greyscale optimization. In the Mg and Na panels some granules 
have larger positive (bright) Dopplershift than others. The re- 
versed granulation sampled in the second row shows the same 



updraft pattern as the granulation in the first row, but the down- 
draft lane pattern is less clear than the updraft granule pattern. 

The observed Mg and Na Dopplergrams in the third row have 
essentially three shades. The bright and dark patches (apart from 
the magnetic concentrations which appear as small black grains) 
correspond roughly to the areas with bright and less bright gran- 
ules in the upper rows, respectively. The intermediate dark-gray 
background consists of unsharp granules (more clearly when 
we flip rapidly through adjacent Dopplergrams on our screen). 
These are not seen in the Ca core Dopplergram, in which the 
patches also have different patterning as well as finer texture. In 
this Dopplergram black grains appear also away from the mag- 
netic concentrations. The strong updraft in the granule near the 
center of the field is visible as a roundish bright patch in the Mg 
and Na cores but not in the Ca core. 

The simulation Dopplergrams also show the granulation pat- 
tern as dominant in the first row, a much fuzzier granulation pat- 
tern with the lanes washed out in the second row, and larger- 
scale patch patterns in the third row. The fine texture in the latter 
is similar for the Mg and Na lines, less so for the Ca line. The 
dark gray shade again corresponds to underlying granules in Mg 
and Na but not in Ca. They are less distinct than in the Mg-core 
brightness images in Figure [6] 

When we scan through all observed Dopplergrams, bright 
updraft patches appear first at AA as 8, 6, 15 km/s for the Mg, 
Na, Ca lines, respectively. Towards the line centers they become 
more opaque. In the simulations they appear only at AA « 4, 5, 
5 km/s, respectively. The simulation Ca Dopplergrams show a 
very flat gray, nearly granule-free low-Dopplershift intermediate 
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Fig. 9. Sample profiles of the cores of the Mg line (top row), Na line (second row), and Ca line (third row), for every 50th pixel 
in x and y in Figures [4] and [5] respectively. First column: observed profiles. Second column: original synthetic profiles from the 
simulation. Third column: synthetic profiles spectrally smeared with the CRISP transmission profile. Fourth column: additional 
smearing with the estimated SST point spread function. The smearing was applied to the computed intensities, not the brightness 
temperatures. Final column: mean profiles. The solid curves are the observed full-field means, the dotted curves the simulation 
means without smearing and with CRISP and SST smearing (the latter does not affect the spatial mean). The dashed curves are the 
FTS atlas profiles without and with CRISP smearing. All panels cover the same coordinate ranges. 



scene over AA — 1 1-6 km/s that is not sampled in the observa- 
tions. 

The more numerous magnetic concentrations in the simula- 
tion also stand out as black blobs in the lower Mg and Na panels. 
The dark patches in the Ca Dopplergram in the last panel have 
some but less good correspondence to the strongest fields. 

4.4. Core asymmetries 

The scatter diagrams in Figure [8] demonstrate the field- 
Dopplershift correlations in more detail. They plot the core 
asymmetry ratio V/R = t^^/T^ against field strength per 
pixel, both for the observations and for the simulation. (We use 
symbo l V for "violet" to corr espond to the Ca n H & K practice 
of e.g. JCram & Darnel (19831) and to avoid confusion with the B 
for magnetic field). The sampling wavelengths for Mg and Na in 
the first and bottom rows are those that show the tightest high- 
high correlations (line-up of the two moment curves). These lie 
further out in the observed profiles than in the synthetic profiles. 
The observed Mg and Na diagrams in the first row show signifi- 
cant correlation between large V/R and large field strength. The 
observed Ca line shows similar correlation at double AA from 



line center. The diagrams in the center row show that the high- 
high correlations remain similar close to line center (tilt of the 
horizontal first-moment curve), but many more weak-field pix- 
els now contribute larger V/R. For all three lines the low-field 
mountain spreads out in V/R, with an extended distribution tail 
for the Na and Ca lines. 

The scatter diagrams from the simulation (bottom row of 
Figure [8} have very different field distributions corresponding 
to Figure [3] They show similar correlations as the observations 
in the top row, but only close to line center. The V/R spreading 
and distribution tails are more like the observed ones near line 
center (center row). 

4.5. Line profiles 

Figure [9] concludes our parallel inspection of observation and 
simulation results. It shows random samples of observed line 
profiles in the first column, random samples of the simulated 
line profiles in the next three columns, and spatially averaged ob- 
served and simulated profiles together with FTS-atlas profiles in 
the last column. The spread between the sample profiles is much 
larger for the simulation than for the observation, but becomes 
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Fig. 10. Strong lines in the red part of the solar spectrum. Disk- 
center profiles of the Na line (Nai 5895.940 A; dotted), the Mg 
line (Mg i 5 172.698 A), the Ca line (Can 8542.144A) and Ha 
(Hi 6562.808 A, dotted) from the F TS atlas are shown on com- 
mon scales. The atlas calibration by Neckel & Labs ( 1984) was 
used to convert the intensities into formal brightness tempera- 
tures. The wavelengths were converted into formal Dopplershifts 
from line center, with redshift positive. 
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Fig. 11. Disk-center profiles of the five strongest Ca n lines from 
the FTS atlas. Format as for Figure [TOl 



similar with the CRISP and SST smearing. The simulated cores 
rise considerably with the CRISP smearing but do not become as 
shallow as the observed ones. The largest discrepancy is obvious 
in all simulation panels, and also in the last column: even after 
the CRISP smearing, the simulated line cores are narrower than 
the observed ones, not much for the Na line but markedly for the 
Ca line. This difference matches the differences in core-structure 
onset with AA and core asymmetry sampling noted above. 

5. Interpretation 

5.1. Solar atlas profiles 

We start with a comparison of the spatially-averaged profiles 
of the three lines. Figure [10] combines them with Ha in a 
common-scale display of their profiles in the FTS atlas. The Na 
line is the deepest, so simple LTE modeling with a radiative- 
equilibrium model (as remains common practice in stellar spec- 
troscopy) would declare this the highest-reaching one, sampling 
the coolest upper-atmosphere layers. However, at the same time 
it is the weakest in terms of equivalent width. The tabulation of 
iMoore ei" al. (1966) gives only 564 mA for the Na line, com- 
pared to 1259 mA for the Mg line and 3670 mA for the Ca line. 
The spread of damping wings also suggests that the Mg line is 
stronger than the Na line, even though its core reaches less deep. 
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Fig. 12. Solid: the three lines in the FTS atlas, as in Figure [TOl 
Dotted: emergent line profiles for the FALC model in the same 
units. 

Ha (equivalent width 4020 mA) is shown for comparison. Its 
core is much wider than the others and even shallower than the 
Mg line. The large core width comes from large thermal broad- 
ening at th e small hydrogen m ass, making it a good temperature 
diagnostic (ICauzzi et al.l 20091). 

The cores of the Na and Mg lines appear similar apart from 
the larger Na line depth and the deeper onset of damping wings 
in the Mg line. The Ca core is markedly wider and displays 
marked asymmetry, whereas the other lines show good sym- 
metry. That this is not due to blending is shown by Figure QT| 
which combines the five strongest solar Ca n lines in the same 
comparison-enabling format. All five display similar red-blue 
core asymmetry. For H & K it appears as the Kxv/Kir asymmetry 
between the slight emission features that set H & K apart from all 
other Fraunhofer lines in the visible as being the only ones not 
showing a simple absorption bell shape. The H & K asymme- 
tries in spatially-averaged profiles a re mostly due to H2V and K2 V 
internetwork grains (cf. review by iRutten & Uite nbroek 1991 ). 
These have been explained by Carlsson & Stein (119941 [19971) 
as acoustic shock interference, taking place in quiet-Sun clapo- 
ti sphe ric gas below the fibrilar Ha chromosphere (cf. iRuttenl 
1995 ). Si milar but field-guided shocks in near-networ k dynamic 
fibrils (lHansteen et al.l 20061; iDe Pontieuet al. 2007) may con- 
tribute similar asymmetry to the atlas profiles. Figure HTIsug- 
gests that the Can infrared triplet lines gain their asymmetry 
likewise. We return to this issue in Section 15.51 and conclude 
from this section that: 

- the Mg and Na lines have similar cores but different wings; 

- the imaging asymmetries in Figure Q] are not accompanied 
by obvious Mg or Na mean-profile asymmetries; 

- only the Ca line displays marked mean-profile asymmetry, 
which it shares with the other strong Ca 11 lines and therefore 
may be due to clapotispheric shocks above the photosphere. 



5.2. Formation in the FALC atmosphere 

After LTE Eddington-Barbier interpretation, the next stage in 
solar spectrum modeling is to relax the assumption of LTE 
but maintain the assumption of hydrostatic equilibrium in a 
plane-parallel simplification. Figure Q~2]shows results from such 
standard NLTE mode ling using the standard FALC model of 
iFontenla et al.l (19931) together with our model atoms in the 
MULTI spectral synthesis code. The FTS atlas profiles of the 
three lines are added for comparison. The pertinent photospheric 
and low-chromosphere parts of FALC are ide ntical to those 
of the quiet-Sun models of [Avrett (1985) and iMaltbv et al.l 
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Fig. 13. Formation paramete rs of the three lines in the stan- 
dard ID model atmosphere of iFontenla et al.l ( 1993) against ge- 
ometrical height. The symbols mark the heights where the total 
optical depth at line center reaches 0.3, 1 and 3, respectively. 
Top panel: NLTE population departure coefficients b for the Mg 
line {dashed), Na line {dotted) and Ca line {solid). The optical 
depth marks are superimposed on the lower-level curves. Middle 
panel: total (opacity-weighted line plus continuum) source func- 
tion in the form of a formal temperature r exc {dashed) and mean 
radiation field in the form of a formal temperature T m ^ {solid) 
for each of the three lines. The gas temperature r gas is added for 
comparison. Bottom panel: scaled height variation of the total 
(line plus continuum) opacity per gram k at the center of each 
line. The curve coding is the same as in the top panel. Note the 
compressed scale for the Ca line; its peak at z-Q is 24 times the 
corresponding peak for the Mg line. 



( 198 6) and represent an u pdate of the canonical VAL3C model 
of Vernazza et al. I dI95H They were defined to fit solar con- 
tinua, not these particular lines, but Figure [l2]shows that they do 
a rather good job for these as well. Our NLTE FALC modeling 
therefore serves well to discuss basic formation aspects of our 
lines. 

Figure[13]diagnoses the FALC formation of the three lines by 
plotting various line formation parameters against height. The 
first panel plots population departure coefficients b = «/«lte 
where n is the computed population and «lte the population that 
would follow from the Saha and Boltzmann equations from the 



abundance, gas temperature and electron density. They demon- 
strate that LTE holds to within 20% in the first few hundred kilo- 
meters of the photosphere. Above the FALC temperature min- 
imum (shown by the r gas curve in the center panel) the lower- 
level curves show steep increases with height that mark NLTE 
underionization in pertinent ionization edges. In these the ion- 
ization is dominated by photospheric radiation from below with- 
out following the FALC chromospheric temperature rise. Ca n 
photo-ionization produces Ca lower-level overpopulation only 
at larger height. 

Deeper down radiative overionization causes small lower- 
level Mg i and Na i underpopulation dips in the upper photo- 
sphereQ. They are less deep than most at omic species had in the 
VAL3C model (e.g., Fe i in iRuttenl 19881) for two reasons. First, 
the FALC upper photosphere has a less steep temperature de- 
cline, and second, for th ese levels such dip s are filled in by pho- 
ton suction, as shown bv lBruls et al . ( 1992) forNa i. This process 
of line-driven recombination also o perates in Mg i and ca uses the 
Mgi 12-micron emission features (ICarlsson et al.l 1992|l. It pro- 
duces the early onsets of the overpopulation rises in Figure [13] 
earliest for the Na line. 

For each line the upper-level departure coefficient (lower 
curve of each pair in the top panel) drops significantly below the 
lower-level one through resonance scattering. The dominance of 
scattering is obvious from the center panel in which the three 
source functions are close to the angle-averaged radiation field 
(in r exc and T a d incarnations to remove Planck function sensi- 
tivity differences). Since S l ~ {b t ,/bi)B (in the Wien limit), the 
divergences between the lower- and upper-level curves in the top 
panel define the divergences between the r g as and r exc curves in 
the center panel. lUitenbroek & Brulsl t 1992) have shown that the 
textbook case of two-level scattering is an excellent approxima- 
tion for the Na line which is a resonance line. The others are not 
and may have larger multi-l evel coupling to the temperature, for 
the C a line via H & K (see lOwocki & Auerl 1980; lUitenbroekl 
19891). 

The Na source function drops deepest because it uncouples 
earliest from the Planck function and feels the FALC chromo- 
spheric temperature rise least. The Ca source function retains a 
hump from this rise and would display small inner-wing emis- 
sion features in FALC modeling if the latter wouldn't apply the 
considerable VAL3C micro turbulence. 

The NLTE source functions and r scaling explain the com- 
puted FALC line profiles in Figure Q~2] Application of the 
Eddington-Barbier approximation to the t- 1 marks in the cen- 
ter panel recovers the computed line depths. 

The bottom panel of Figure [13] shows the buildup of the t 
scales, which differs between the three lines (note the scale fac- 
tors needed to overplot these curves within one panel). Plotting 
the opacity per gram takes out the density dependence. The curve 
patterns are the product of the Saha and Boltzmann population 
formulas, the shape of the line extinction profile, and the lower- 
level population departures bi in the top panel. 

Let us first discuss the k curve for the Na line (dotted). The 
lower-level population/density ratio ni/p of the Na line in LTE 
is set by the ionization equilibrium since this line is from the 



1 More recently. IFontenla ~et al. (2009) claimed that the much deeper 
upper-photosphere b dips for Si n and the other electron-donor species 
in their newer models are due to irradiation from above, not from be- 
low, but this is a misinterpretation due to using the M enzel definition for 
depart ure coefficients, see page 3 7 of Rutten (2003J). We use the defini- 
tion of Wijbenga & Zwaan ( 1972) which normalizes b by the element 
abundance instead of the next-ion density. 
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ground state and has no Boltzmann excitation sensitivity to tem- 
perature. This ratio (not shown) has an initial rise due to de- 
creasing ionization in the initially steeply declining temperature, 
followed by a long decline due to larger ionization at lower elec- 
tron density. The photospheric wiggle pattern in the actual Na- 
line k curve is added by the b\ variations in the top panel and 
by changes in the line profile shape. The steep b\ rise halts the 
ionization decay near z = 0.7 Mm. 

The Mg line has three times larger opacity in the lower pho- 
tosphere but a much larger drop across the upper photosphere 
due to the sensitivity of its Boltzmann excitation to the declining 
temperature. The Na and Mg FALC opacities are nearly identical 
above /i = 0.5 Mm. 

The Ca line is from the dominant calcium ionization stage; 
ionization of Ca n to Ca m affects At a from z = 0.9 Mm. In the 
photosphere this line has no ionization sensitivity, only excita- 
tion sensitivity to the declining temperature. This sensitivity is 
smaller than for the Mg line (1.7 eV versus 2.7 eV excitation 
energy). 

The wiggles in the deepest layers are set by the steep out- 
ward decrease of the electron/gas density ratio due to hydrogen 
recombination. It affects the ionization and also the balance be- 
tween Stark and Van der Waals broadening in setting the broad 
wings of the Ca and Mg lines. The damping constant of the Mg 
line is about 1 .7 times that of the Na line across the FALC pho- 
tosphere and much larger in the Stark domain below h = Mm. 

We conclude from our FALC modeling: 

- the Mg line has much larger opacity and somewhat larger 
damping than the Na line in the lower photosphere, hence its 
more prominent damping wings in Figure [10] and its twice- 
larger equivalent width; 

- its Boltzmann sensitivity to the declining temperature across 
the photosphere reduces the Mg line opacity to that of the Na 
line, so that they have nearly identical t— 1 heights in FALC; 

- although the FALC chromosphere is not a realistic model of 
the actual chromosphere, this higher-layer opacity equality 
suggests that chromospheric fibrils can have similar opaque- 
ness in the Mg and Na lines, as observed in Figure Q] 

- the various contributions to the k variations around /; = 65 km 
in Figure[T3]with their different temperature and density sen- 
sitivities produce the contrast differences along the first and 
last rows of Figure [6] For example, the Ca line has lower 
contrast because it is not sensitive to ionization while its 
B oltzmann excitation cau ses contrast reduction (cf. Figure 7 
of lLeenaarts et al.l20 06b ). The Mg and Na lines gain contrast 
from the k sensitivities to ionization that produce dips at this 
height in Figure [T3l 

- the Mg-line opacity has larger temperature sensitivity than 
the Na-line opacity due to its excitation energy. This may 
contribute to the larger clapotisphere transparency for Mg 
that is seen in the center row of Figure |5] 

- the Ca line has much larger opacity across all of FALC, 
which explains that it also maps chromospheric fibrils bet- 
ter than the other two lines. In the upper photosphere the Ca 
wings sample identical scenes as the other two lines but fur- 
ther out from line center, as seen in Figures [5]-|6] 

- difference in resonance scattering explains that the Na line- 
center image in Figure |4] appears the fuzziest, the Ca line- 
center image the sharpest. Similar sharpness difference is 
seen between the line-center Dopplergrams in Figure [7] also 
for the simulation. 
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Fig. 14. Representative observed Na-line profiles to illustrate 
the cause of imaging asymmetries in granulation and reversed 
granulation. Solid: pixel at the center of a granule. Dashed: pixel 
in an intergranular lane. Dotted: spatial average. 



5.3. Granulation and reversed granulation 

Comparison of the second and next-to-last rows in Figures |4|- 
[6] confirms the imaging asymmetries in Figure Q] Figure [TO] 
shows that it does not have to do with mean-profile asymme- 
try. Figure [13] demonstrates that the Na line and Mg line have 
near-equal FALC core formation. Therefore, the reason for the 
imaging asymmetry must be dynamics that is not included in 
the FALC modeling. The Dopplergrams in Figure ^demonstrate 
that the reversed-granulation scenes in Figures HJE|are fully set 
by granular dynamics. The latter causes the imaging asymme- 
tries. 

Figure [14] is a simple demonstration for the Na line follow- 
ing th e cartoon for Ban 4554 A in Figure 6 of ISiitterlin et al.l 
(2001) with act ual data. (Sim ilar cartoons were already shown 
in Figure 4 of lEvansI 19641) The lane pixel has a redshifted 
profile that is darker than the mean profile, especially in the 
outer wings. The granule pixel has a profile close to the spa- 
tial mean but slightly blueshifted and slightly brighter in the 
outer wings. A filtergram near AA = -14 km/s as the first Na 
panel in Figures[4}]6]shows "flattened" granulation with smaller 
contrast than between nonshifted profiles. The reverse holds for 
the opposite sampling at AA = +14 km/s. This explains that in 
Figures [4}]6] the top Na panel shows lower granulation contrast 
than the bottom Na panel. The corresponding Dopplershift pan- 
els in Figure [7] show the granulation even clearer, both for the 
observation and the simulation. 

The same comparison for the Ca line shows no such differ- 
ence, the corresponding observed Dopplergram only noise. The 
reason is simply that the comparable scene (or height) sampling 
takes place much further out in the extended damping wings of 
the Ca line, where these have a shallow slope so that profile shift 
has less effect. This is also the case for the Mg line which differs 
from the Na line in possessing extended wings (Figure [TOt. 

Let us now make such comparisons closer to the core. The 
blue wings of the lane and granule profiles in Figure [14] cross, 
giving granulation contrast reversal which does not occur on the 
red side. These two pixels show reversed granulation with re- 
versed imaging asymmetry. The effect is again larger for steeper 
line flanks. The corresponding sampling temperature of about 
5000 K lies above the onset of the Mg and Ca wings in Figure [T2l 
Therefore, the striking imaging asymmetry between the wings 
of the Na line in Figure Q] is also simply due to their steepness 
(Figure [TOt . The same asymmetry should appear in the inner 
core of other narrow but somewhat weaker lines; indeed it is very 
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pronou nced for Fei 7090.4 A in Figure 5 of Ijanssen & Cauzzil 

MR 

The diffuse-gray lane patterns in the middle-row 
Dopplergrams in Figure UJ correspond well to the vertical 
velocity pattern at z = 0.3 Mm in the last row of Figure 5 of 
iLeenaarts et al.l (20101) and qualitative l y to th e similar pattern at 
t = 0.1 in Figure 3 of ICheung et al.l (2007) Such simulations 
also show slight temperature contrast reversal between granules 
and lanes at such heights, in our case illustrated for z - 0.3 Mm 
in the second row of Figure 5 o f ILeenaarts eTaD (20101) It was 
already explained by iNordlundl ( 1984bl) on pages 15-18 in his 
informative analysis of his pioneering granulation simulation 
as due to radiative reh eating along bend-ov er flow trajectories, 
as those in Figure 6 of ICheung et al. (2007) The Dopplershift- 
brightening pattern combination produces slender near-lane 
brightness features as seen in the second row of Figure |4] They 
correlate only roughly with underlying intergranular lanes, best 
at a delay of 1-2 minutes. In addition, at larger height gravity 
wave interference takes over fro m direct granular o vershoo t and 
reduces the correspondence ( cf. iRutten et al . 2004 L eenaarts & 
Wede meyer-Bohm 2005 ; Jan ssen & Cauzzil 20061, IStraus et all 
20081). 

The conclusion from this section is that the imaging asym- 
metries between the Na wings and between the Mg and Na lines 
in Figure[T]arise from the correlated brightness and Dopplershift 
patterns in the solar granulation and its overshoot, to which the 
Na line is more sensitive because it does not have strong damp- 
ing wings. 

5.4. Line core scenes 

All three lines are scattering lines, the Na line the worst. Their 
core intensities are set by radiation escape in deeper layers than 
the last scatterings. The latter set the Dopplershift encoding (as- 
suming com plete redistribution which i s in order even for the 
Na Line, see lUitenbroek & Brulsl 19921). This sampling differ- 
ence is most clearly demonstrated by the simulated Mg core, 
which shows regular granulation much clearer in intensity than 
in its Dopplergram. The difference also explains that the large 
image variations between the three central rows of Figure @J[6] 
are largely produced by the Dopplershifts in the bottom rows 
of Figure [7J Therefore, taking Dopplergrams represents a much 
better use of these lines than taking intensity filtergrams. 

The rough correspondence of the bright and dark Mg and Na 
Dopplergram patches in the third row of Figure UJ with the ar- 
eas of brighter and less bright granules in the first row suggests 
that these are p-mode oscillation patches that result from global 
interference patterning, with larger D in the blueshifted phase. 
The lambdame ter phase difference measurements for Na i Dt of 
iDeubner et al.l ( 19961) indeed show oscillation coherence be- 
tween the outer wing and core, but only for the five-minute 
component (last panel of their Figure 6). For increasing height 
the oscillation dominance shifts from five-minute evanescent 
global mode interference to more local three-minute propagating 
waves, as seen already in the somewhat deepe r layers sampled by 
ultraviolet continua (e.g.. lKriiger e t al. 2001 ), and then to three- 
minute shocks that are seen in the Ca line (next section). The 
lack of coherence at higher frequency in Figure 6 of Deubner 
et al. (19961) suggests that the Mg and Na core Dopplergrams 
already have some sensitivity to these shocks. 

The simulation does not contain global p-mode oscillation 
interference patterning but has its own box modes and gener- 
ates its own propagating waves. The fine structure across the 
last row of Figure [7J shows that all three lines Doppler-sense 



the very inhomogeneous shock-ridden clapotisphere that results, 
even though it appears too transparent in the simulate d Mg line 
(Section©. We refer to Leenaarts et al. (20091 20101) for perti- 
nent displays. In particular, Figure 9 of the latter paper contains 
informative formation breakdown diagrams for the Na and Ca 
lines. Comparing these to Figure[13]demonstrates enormous line 
formation variations with respect to FALC modeling. 

The conclusion of this section is that the line cores primarily 
sample wave patterns that likely include shock interference. 

5.5. Magnetic concentrations and shocks 

There are only a few magnetic concentrations in the observed 
subfield, more in the simulation. In Figures |4]-|6] they stand out 
as bright features in the outer line wings, often brighter than 
the brightest granules because lack of collisional damping in the 
partially evacu ated concentrat ions reduces the wing opacity (cf. 
Leenaarts et al. 2006a| 2006bl). 

Nearer the line centers they no longer stand out as bright 
points but rather as profile asymmetry features, producing small, 
roundish, dark features in the Mg and Na Dopplergrams in 
Figure [7] Figure [8] uses the V/ R asymmetry rat io following the 
example for Can Hhy grains of ICram & Dam e (1983) The for- 
mation diagrams of such grains in lCarlsson & Steir3 ( 1997) show 
that JLy grain asymmetry results from line-of-sight integration 
through cool postshock downflows around h — 1.5 Mm and 
grain-producing upcoming hot fresh shocks around h — 1 Mm. 
The same pattern occurs for rarefied magnetic concentrations 
at heights of only a few hundred kilometers, where the outside 
scene is still reversed granulation. The ubiquity of photospheric 
shocks in an d near isolated ma gnetic concentrations was first 
suggested bv lSteiner et al . ( 1998) Their signature in the Na core 
was also noted bv ljess et al.l (120101)1 

A detailed demonstration is given b y the Na-line formatio n 
diagrams in the first quartet of Figure 7 of lLeenaarts et al. 1 (20101) 
which analyze the magnetic concentration at (x,y) = (8.2,7.2) 
in the simulation snapshot. It has a shock at z » 0.3 Mm 
that produces an emission peak just blue of Na line center, in 
the same manner as a H2V grain but much deeper. The post- 
shock downdrafts from shocks that previously passed, now at 
z ~ 0.5 and 1.0 Mm, produce a red-shifted line core which is 
very dark because the strongly scattering Na source function de- 
couples from the Planck function already at h * 0.2 Mm and 
shows no sensitivity to the actual high-temperature plateau over 
h ss 0.3 - 0.7 Mm. Such raised temperature plate aus in magnetic 
concentrations are attributed to Joule heating by ICarlsson et al.l 
(2010) They are similar to the FALC chromospheric tempera- 
ture plateau but occur at photospheric heights. 

Note that when such a spectral emission peak blueshifts 
into an inner-wing sampling passband, the Dopplergram sign of 
D = {T l b ed - T™ ue )/(r b ed + r£ lue ) becomes the same as for a 
redshifted absorption core. K2y-fike grains can therefore appear 
black in our D definition. Similarly, the asymmetry ratio V/R of 
Figure[8]may exceed unity both due to fresh-shock bright updraft 
and to dark post-shock downfall, i.e., for both phases of shock 
modulation and especially when these combine along the line of 
sight. 

The Ca image just blue of line center in the third row of 
Figure|4]does show non-magnetic bright grains that represent the 
Ca-line analog to higher-formed internetwork K2V grains. Their 
visibility is to be expected from the similar core asymmetries 
in all five Ca 11 lines in Figure QT| They are indeed black in the 
corresponding Dopplergram in Figure UJ The center-row Ca di- 
agram in Figure [8] suggests larger spatial filling factor for these 
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Fig. 15. Comparable intensity images from 
Hinode in the Mg line (upper half) and Na 
line (lower half). Each panel is independently 
grayscaled for optimum contrast. The subfield 
has the same size as in Figures 0-[7] Along 
rows the images cover the same scene, with 
about 32 s delay between the two lines. The 
viewing angle differs between rows, with 
p = 0.35, 0.50 and 1.00 as specified. The limb 
is to the top. The first two columns are images 
in the blue and red wing at A A = +240 mA 
for the Mg line, A/1 = ± 160 mA for the Mg 
line. The third column has the corresponding 
(^ed - 4iu e )/<7red + /blue) Dopplergrams. 
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shocks. The simulation diagram even has a secondary contour 
peak for them. 

Outside magnetic concentrations the Na and Mg lines are too 
weak to produce internetwork grains. Nevertheless, the Na-core 
V/R distribution tail in the central diagram of Figure |8]suggests 
that the Na line senses internetwork shocks similarly to the Ca 
line. 

The conclusion of this section is that Na and Mg core 
Dopplergrams or asymmetry images turn out to be quasi- 
magnetograms by mapping photospheric shocks in or near 
strong-field concentrations. Outside these, only the Ca line 
shows mean-profile asymmetry from clapotispheric internet- 
work shocks. 

6. Center-to-limb variation 

Figure Q3] presents a selection from the Hinode observations, 
added to show how limb ward viewing affects the three lines. 
The AA values are selected for scene similarity at pi = 1 . These 
images (bottom rows of each half) have lower quality than the 
SST/CRISP images but again show similar imaging asymmetry 
as in Figure [T] similar reversed granulation in the blue Mg and 
Na wings and the red Mg wing, granular-like patterning instead 
in the red Na wing. The p-l Dopplergrams also show granula- 
tion, with larger contrast in the Na line. 

At p = 0.5 the pattern repeats but with less discrepancy for 
the Na red wing, less difference between the Dopplergrams, and 
general foreshortening. The decrease in imaging asymmetry cor- 
responds to smaller Dopplershifts from the vertical motions. 

Closer to the limb a different pattern of horizontal stripings 
appears, similarly in all p = 0.35 panels and clearest in the Mg 
Dopplergram. It looks rathe r like the supergranula tion in the dis- 
covery plate in Figure 6 of iLeighton et al. (1962) but of course 



the scale is much smaller: these are horizontal outflows in gran- 
ulation cells rather than in supergranulation cells. At this for- 
mation height they probably correspond to the bend-over flow 
lines of lCheung et al.1 (2007) In the Na Dopplergram some dark 
patches are visible which are likely clapotispheric clouds that are 
not visible in the Mg line due to lack of Boltzmann excitation. 



7. Discussion 

In Figures [4}j7] we chose to jump from deep outer-wing forma- 
tion (h as 150 km) in the second rows to line-core formation 
in the next rows to avoid figure clutter and to skip observation- 
simulation sampling differences set by the too narrow simu- 
lated Ca and Mg cores. In the observations the intervening 
scenes, best seen per Dopplergram, are simply the combination 
of the (reversed) granulation and overlying larger-scale oscilla- 
tion patches. The cloud-like appearance of the latter suggests 
bimodal formation, with the deeper layers well described by 
Eddington-Barbier sampling near t = 1, the upper layers better 
described by Schuster-Schwarschild clouds with Dopplershift as 
the most i mportant cloud parame ter. The inner-core panels in 
Figure 6 of iDeubner et al.l (19961) confirm the notion of cloud- 
like common Dopplershift at all frequencies. In the simula- 
tion the visibility of a gray low-Dopplershift intermediate back- 
ground scene that appears especially in the outer Ca core is likely 
due to incomplete masking by clapotispheric clouds that are too 

transparent. 

Cauzzi et al. (20091 20081) suggested such bimodal formation 
for chromospheric fibrils in the Ca line. These indeed blanket 
the inner Ca core even more when there is more activity than in 
the very quiet area analysed here. In the B aimer lines they eve n 
tend to o bscure the reversed granula tion (tLeenaarts et al.l 2 006b ). 
Note that lEvans & Catalanol (1972!) already reported not finding 
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Fig. 16. Dopplergrams from the simulation as in the bottom row of Figure |7]but closer to line center and with only CRISP spectral 
smearing, no spatial PSF smearing. 



negative "oddities" for H/3 (but note also that lThomasI ( 1972|) as 
referee blamed fibril source functions, not opacities). 

The simulation snapshot and line synthesis do a good job in 
reproducing and so helping to explain the imaging asymmetry 
between the Na line and the other two lines, but they do not re- 
produce the observations in the line cores very well. These are 
not the main topic here so we discuss th ese discrepan cies only 
briefly, again referring to Leenaarts et al. (2009 2010 1 for more 
detail. The major ones are the larger visibility of granulation at 
the center of the Mg line, the narrower core width of the Ca line 
in the simulation, and the differences in core asymmetry sam- 
pling in Figure [8] 

The clapotispheric cool-cloud transparency in the simulated 
Mg line is likely an artifact due to underestimation of the Mg i 
density, for the same reason as the und erestimation of the Na i 
density in these cool clouds discussed in lLeenaarts et al.l (20101): 
since in the simulation the ionization equilibria are assumed to 
be in LTE, the electron-donor elements (Si, Fe, Al, Mg) are less 
ionized than they would be if the actual domination by photoion- 
ization was properly accounted for. The resulting underestima- 
tion of the electron density is 1-2 orders of magnitude, higher 
up even 3-4 if one compares with time-d ependent hydrogen ion - 
ization balancing in cool post-shock gas dLeenaarts et al.l 20071). 
In contrast, the synthesis codes implement photoionization prop- 
erly, so these yield higher ionization rates, but the recombination 
lacks electrons so that the balance shifts to the ion state and the 
line opacity is severely underestimated. Thus, the cool clouds are 
too transparent in the simulated Na and Mg lines, more so in the 
latter due to its 2.7 eV excitation energy. Possibly the Ca line is 
similarly affected higher up in cool clouds by similarly overes- 
timated ionization from Ca n to Ca m. Such transparencies may 
contribute to the differences in Figure [8] 

The to o narrow core wi dth of the Ca line was noted and dis- 
cussed by iLeenaarts et a l. (2009) The apparent need for "mi- 
croturbulent broadening" in higher layers suggests a lack of 
numerical resolution, just as the solar granulation needed mi- 
croturbulent fu dging until sim ulation resolution did away with 
that (Nordlund 19801 1984b). In deed, a t est w ith the higher- 
resolution simulation of ICarlsson et a l. (2010) shows larger 
broadening. The assumption of instantaneous LTE ionization 
balancing not only causes underestimation of the electron den- 
sities in cool clouds and corresponding underestimation of the 
opacitie s of our three li nes, but also of the clapotispheric dy- 
namics dLeenaarts et al.l 20071). 

These deficiencies suggest that the actual clapotispheric con- 
tributions to these line cores may be even more wildly struc- 
tured than Figures [6}]7] suggest already. Figure [16] displays 
innermost-core Dopplergrams from the simulation after only 
spectral smearing with the CRISP transmission functions. The 



spectral smearing affects the core depth considerably (Figure |9) 
but the Dopplergram scene keeps its structure. The scene is a vi- 
olent one with sharp bright fronts, dark cool clouds, and black 
magnetic concentrations. Comparison of Figure [T6lto the Na 
line-ce nter Dopplershift panel in Figure 7 of ILeenaarts et al.l 
(20101) shows that the CRISP-smeared Dopplergram gives a 
good rendering of the actual Dopplershifts, with more empha- 
sis on magnetic concentrations. Comparison of Figure Q~6] to 
the center rows of Figure |5] shows that the simulated line cores 
are so much dominated by Doppler modulation that inner-core 
Doppler measurement is mandatory in using these lines as upper- 
atmosphere diagnostics. Comparison of Figure [16] to the line- 
core rows of Figure Qsuggests that higher spatial resolution than 
in our observations is highly desirable. 

A corollary is that searches for global g-modes and so-called 
"chromosphere seismology" using full-disk resonance-cell sam- 
pling of the Na line as wit h GOLF-NG (e.g.. iTurck-Ch ieze et al. 
2006 JSalabert et al. 2009 ) will suffer noise from reversed granu- 
lation just as classical helioseismology suffers noise from gran- 
ulation. Magnetic concentrations contribute much noise (liter- 
ally) by their shocks. GOLF-NG's 15 passbands, spread over 
AA = +9 km/s, have FWHM = 30 mA, half CRISP's pass- 
band in this line and clearly narrow enough. Interpretation of 
such multi-passband oscillation sampling that relies on one- 
dimensional height-of-formation interpretation, for example to 
diagnose upward propagation from inward phase difference, is 
likely to fail since the clapotispheric signals are better described 
as clouds of varying opacity at varying height of which the vary- 
ing Dopplershifts act as shutters obscuring the inner line wings. 



8. Conclusion 

There is a striking imaging asymmetry between Nai Di and 
Mgib2 in Figure [TJ the Mg line shows reversed granulation 
in both wings, but the Na line shows similar reversed granula- 
tion only in its blue wing, normal granulation in its red wing. 
This is due to different sensitivity to the correlated Dopplershift 
and brightness pattern of the reversed granulation in the mid- 
photosphere. The Mg line has larger opacity and damping in 
deeper layers, hence wider damping wings, hence less steep pro- 
file flanks, hence smaller Doppler sensitivity. Both lines sample 
the same layers in both their wings; only the appearance of the 
granulation and the reversed granulation differ between these, 
not their physical nature or sampling height. 

The cores of the Mg and Na lines are formed remarkably 
similarly. Dopplergrams in these represent their primary diag- 
nostic value and can also serve to detect and study magnetic 
concentrations marked by photospheric shocks. 
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Center-to-limb observations with at least the quality of the 
disk-center SST/CRJSP data used here will be well-suited to 
constrain the reversed-granulation domain further and to enlarge 
the testing ground for numerical simulations. 

When the lack of line-core opacity and broadening in the 
simulation is resolved, the obvious next step is to compare time 
sequences of such imaging spectroscopy an d simulations for 
these l ines in Fourier analyses as those of e.g. jFleck & Deubnerl 
( 19891) and lDeubneret al.l (19961) 

The small structural scales that are predicted by the simula- 
tions to dominate at clapotispheric heights define the need for 
high spatial resolution, considerably higher even than reached 
here, as a principal quest in probing this domain of the solar at- 
mosphere. 

Acknowledgements. We thank T.D. Tarbell and P. Siitterlin for the discussions 
that inspired this paper and T. Leifsen, H. Skogsrud and A. Ortiz for obtaining the 
CRISP scan in Figurefjthat initiated these discussions. The SST is operated by 
The Institute for Solar Physics of the Royal Swedish Academy of Sciences at the 
Spanish Observatorio del Roque de los Muchachos of the Instituto de Astroffsica 
de Canarias. The Japanese Hinode mission was developed and launched by 
ISAS/JAXA, with NAOJ as domestic partner and NASA and STFC as inter- 
national partners, and is operated by these agencies in co-operation with ESA 
and NSC. RJR's travel to the SST and to Oslo was funded by the EC through the 
SOLAIRE Network (MTRN-CT-2006-035484). Discussions at NAOJ (Mitaka, 
Japan) while RJR was a visitor there also contributed to this paper, and he com- 
pleted it as a visitor at LMSAL (Palo Alto, USA). As always, our work was much 
facilitated by NASA's ADS. 

References 

Altrock, R. C. & Musman, S. 1976, ApJ, 203, 533 

Avrett, E. H. 1985, in Chromospheric Diagnostics and Modelling, ed. 

B. W. Lites, 67 
Bard, S. & Carlsson. M. 2008, ApJ, 682, 1376 
Botnen, A. 1997, Master's thesis, Inst. Theor. Astrophys. Oslo 
Bruls, J. H. M. J., Rutten, R. J., & Shchukina, N. G. 1992, A&A, 265, 237 
Canfield, R. C. & Mehltretter, J. P. 1973, Sol. Phys., 33, 33 
Carlsson, M. 1986, Uppsala Astron. Obs. Reports, 33 
Carlsson, M., Hansteen, V. H., & Gudiksen, B. V. 2010, MemSAI, 81, 582 
Carlsson, M., Rutten, R. J., & Shchukina, N. G. 1992, A&A, 253, 567 
Carlsson, M. & Stein, R. F. 1994, in Chromospheric Dynamics, ed. M. Carlsson, 

47 

Carlsson, M. & Stein, R. F. 1997, ApJ, 481, 500 

Cauzzi, G, Reardon, K., Rutten, R. J., Tritschler, A., & Uitenbroek, H. 2009, 
A&A, 503, 577 

Cauzzi, G, Reardon, K. P., Uitenbroek, H„ et al. 2008, A&A, 480, 515 

Cheung, M. C. M.. Schussler, M., & Moreno-Insertis, F. 2007, A&A, 461, 1163 

Collados, M. & Vazquez, M. 1987, A&A, 180, 223 

Collet, R., Asplund, M., & Thevenin, F. 2005, A&A, 442, 643 

Cram, L. E. & Dame, L. 1983, ApJ, 272, 355 

De Pontieu, B., Hansteen, V. H., Rouppe van der Voort, L., van Noort, M., & 

Carlsson, M. 2007, ApJ, 655, 624 
Deubner, F, Waldschik, T., & Steffens, S. 1996, A&A, 307, 936 
Deubner, F. L. & Mattig, W. 1975, A&A, 45, 167 
Evans, J. W. 1964, Astrophysica Norvegica, 9, 33 
Evans, J. W. & Catalano, C. P. 1972, Sol. Phys., 27, 299 
Fleck, B. & Deubner, F. 1989, A&A, 224, 245 
Fleck, B., Straus, T, Carlsson, M., et al. 2010, MemSAI, 81, 777 
Fontenla, J. M., Avrett, E. H., & Loeser, R. 1993, ApJ, 406, 319 
Fontenla, J. M., Curdt, W., Haberreiter, M., Harder, J., & Tian, H. 2009, ApJ, 

707, 482 

Hansteen, V. H., Carlsson, M., & Gudiksen, B. 2007, in ASP Conf. Ser. 368, 
The Physics of Chromospheric Plasmas, ed. P. Heinzel, I. Dorotovic, & 
R. J. Rutten, 107 

Hansteen, V. H., De Pontieu, B., Rouppe van der Voort, L., van Noort, M., & 

Carlsson, M. 2006, ApJ, 647, L73 
Janssen, K. & Cauzzi, G. 2006, A&A, 450, 365 

Jess, D. B., Mathioudakis, M., Christian, D. J., Crockett, P. J., & Keenan, F. P. 

2010, ApJ, 719, L134 
Kosugi, T, Matsuzaki, K., Sakao, T, et al. 2007, Sol. Phys., 243, 3 
Krijger, J. M., Rutten, R. J., Lites, B. W., et al. 2001, A&A, 379, 1052 
Leenaarts, J. & Carlsson, M. 2009, in ASP Conf. Ser. 415, The Second Hinode 

Science Meeting: Beyond Discovery-Toward Understanding, ed. B. Lites, 

M. Cheung, T. Magara, J. Mariska, & K. Reeves, 87 



Leenaarts, J., Carlsson, M., Hansteen, V., & Rouppe van der Voort, L. 2009, ApJ, 
694, L128 

Leenaarts, J., Carlsson, M., Hansteen, V, & Rutten, R. J. 2007, A&A, 473, 625 
Leenaarts, J., Rutten, R. J., Carlsson, M., & Uitenbroek, H. 2006a, A&A, 452, 
L15 

Leenaarts, J., Rutten, R. J., Reardon, K., Carlsson, M., & Hansteen, V. 2010, 
ApJ, 709, 1362 

Leenaarts, J., Rutten, R. J., Siitterlin, P., Carlsson, M., & Uitenbroek, H. 2006b, 

A&A, 449, 1209 
Leenaarts, J. & Wedemeyer-Bohm, S. 2005, A&A, 431, 687 
Leighton, R. B., Noyes, R. W., & Simon, G. W. 1962, ApJ, 135, 474 
Levy, M. 1971, A&A, 14, 15 

Lockyer, N. 1868, Procs. Royal Soc. London, 17, 131 

Maltby, P., Avrett, E. H., Carlsson, M., et al. 1986, ApJ, 306, 284 

Moore, C. E., Minnaert, M. G. J., & Houtgast, J. 1966, The solar spectrum 

2935 A to 8770 A 
Neckel, H. 1999, Sol. Phys., 184, 421 
Neckel, H. & Labs, D. 1984, Sol. Phys., 90, 205 

Nordlund, A. 1980, IAU Colloq. 51: Stellar Turbulence, ed. D. F. Gray & 
J. L. Linsky, Lecture Notes in Physics, Springer, 1 14, 213 

Nordlund, A. 1984a, in Small-Scale Dynamical Processes in Quiet Stellar 
Atmospheres, ed. S. L. Keil, 174 

Nordlund, A. 1984b, in Small-Scale Dynamical Processes in Quiet Stellar 
Atmospheres, ed. S. L. Keil, 181 

Owocki, S. P. & Auer, L. H. 1980, ApJ, 241, 448 

Roddier, F. 1981, Prog. Optics, 19, 281 

Rouppe van der Voort, L., Leenaarts, J., De Pontieu, B., Carlsson, M., & Vissers, 

G. 2009, ApJ, 705, 272 
Rutten, R. J. 1988, in IAU Colloq. 94: Physics of Formation of FE II Lines 

Outside LTE, ed. R. Viotti, A. Vittone, & M. Friedjung, ASSL, Springer, 138, 

185 

Rutten, R. J. 1995, in ESA Special Pub. 376, Helioseismology, 151 
Rutten, R. J. 2003, Radiative Transfer in Stellar Atmospheres, Lecture notes 
Utrecht University 

Rutten, R. J. 2010, in Recent Advances in Spectroscopy: Astrophysical, 
Theoretical and Experimental Perspectives, ed. R. K. Chaudhuri, M. V. 
Mekkaden, A. V. Raveendran, & A. S. Narayanan, Astrophys. Space Science 
Procs., Springer, 163, preprint 2009arXiv0912.2206L 

Rutten, R. J., De Wijn, A. G, & Siitterlin, P. 2004, A&A, 416, 333 

Rutten, R. J. & Uitenbroek, H. 1991, Sol. Phys., 134, 15 

Salabert, D., Turck-Chieze, S., Barriere, J. C, et al. 2009, in ASP Conf. Ser. 

416, Solar-Stellar Dynamos as Revealed by Helio- and Asteroseismology: 

GONG 2008/SOHO 21, ed. M. Dikpati, T. Arentoft, I. Gonzalez Hernandez, 

C. Lindsey, & F. Hill, 341 
Scharmer, G. B., Bjelksjo, K., Korhonen, T. K., Lindberg, B., & Petterson, B. 

2003a, SPIE Conf. Series 4853, 341 
Scharmer, G. B., Dettori, P. M., Lofdahl, M. G, & Shand, M. 2003b, SPIE Conf. 

Series 4853, 370 

Scharmer, G. B., Narayan, G, Hillberg, T., et al. 2008, ApJ, 689, L69 
Schussler, M., Shelyag, S., Berdyugina, S., Vijgler, A., & Solanki, S. K. 2003, 
ApJ, 597, L173 

Shelyag, S., Schussler, M., Solanki, S. K., Berdyugina, S. V, & Vogler, A. 2004, 
A&A, 427, 335 

Steiner, O, Grossmann-Doerth, U., Knolker, M., & Schussler, M. 1998, ApJ, 
495, 468 

Straus, T., Fleck, B., Jefferies, S. M., et al. 2008, ApJ, 681, L125 
Suematsu, Y., Tsuneta, S., Ichimoto, K., et al. 2008, Sol. Phys., 249, 197 
Suemoto, Z., Hiei, E„ & Nakagomi, Y. 1987, Sol. Phys., 112, 59 
Suemoto, Z., Hiei, E., & Nakagomi, Y. 1990, Sol. Phys., 127, 11 
Siitterlin, P., Rutten, R. J., & Skomorovsky, V. I. 2001, A&A, 378, 251 
Thomas, R. N. 1972, Sol. Phys., 27, 303 

Tsuneta, S., Ichimoto, K., Katsukawa, Y, et al. 2008, Sol. Phys., 249, 167 
Turck-Chieze, S., Carton, P., Ballot, J., et al. 2006, Adv. Space Res. 38, 1812 
Uitenbroek, H. 1989, A&A, 213, 360 
Uitenbroek, H. & Bruls, J. H. M. J. 1992, A&A, 265, 268 
van Noort, M., Rouppe van der Voort, L., & Lofdahl, M. G. 2005, Sol. Phys., 
228, 191 

van Noort, M. J. & Rouppe van der Voort, L. H. M. 2008, A&A, 489, 429 
Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 
Wedemeyer-Bohm, S. 2008, A&A, 487, 399 

Wedemeyer-Bohm, S. & Rouppe van der Voort, L. 2009, A&A, 503, 225 
Wijbenga, J. W. & Zwaan, C. 1972, Sol. Phys., 23, 265 



